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An exclusive measurement of the excitation function for the dd — > 3 Hep7r~ reaction was performed 
at the Cooler Synchrotron COSY- Julich with the WASA-at-COSY detection system. The data were 
taken during a slow acceleration of the beam from 2.185 GeV/c to 2.400 GeV/c crossing the kinematic 
threshold for the rj meson production in the dd — > 4 Her; reaction at 2.336 GeV/c. The corresponding 
excess energy with respect to the 4 He — rj system varied from -51.4 MeV to 22 MeV. The integrated 
luminosity in the experiment was determined using the dd — > 3 Hen reaction. The shape of the 
excitation function for the dd —¥ 3 Hep7r~ was examined. No signal of the 4 He — rj bound state 
was observed. An upper limit for the cross-section for the bound state formation and decay in the 
process dd — >• ( 4 He — rj)bound ~ > 3 Hep7r~, was determined on the 90 % confidence level and it varies 
from 20 nb to 27 nb for the bound state width ranging from 5 MeV to 35 MeV, respectively. 
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PACS numbers: 21.85.+d, 21.65.Jk, 25.80.-e, 13.75.-n 
Keywords: mesic nuclei, rj-mesic nucleus, r] meson 



INTRODUCTION 



EXPERIMENT 



Neutral mesons such as e.g. 77, K, u, rj' can potentially 
form bound states with atomic nuclei. In this case the 
binding is exclusively due to the strong interaction and 
the bound state or mesic nucleus - can be considered as 
a meson captured in the mean field of the nucleons. Due 
to the strong attractive 77-nucleon interaction P, Hj, the 
?7-mesic nuclei are some of the most promising candidates 
for such states. 

Experimental confirmation of the existence of 77-mesic 
nuclei would be interesting on its own but it would be also 
valuable for investigations of the n — N interaction and for 
the study of in-medium properties of the N* resonance |S] 
and of the 77 meson [4]. It could also help to determine 
the flavour singlet component of the 77 wave function [|| . 

The existence of 77-mesic nuclei was postulated in 1986 
by Haider and Liu [6] . Experimental searches have been 



performed by several past experiments [7 - 11 
ing investigations continue at COSY |12Hl7j . JINR If 
J-PARC 0, MAMI [Hi and are planned at GSI |2jJ . 
Many promising indications where reported, however, so 
far there is no direct experimental confirmation of the 
existence of mesic nuclei. 

In the region of the light nuclei systems such as e.g. 
77-He the observation of a strong enhancement in the to- 
tal production cross-section and the phase variation of 
the scattering amplitude in the close-to-threshold region 
provided strong evidence for the existence of a pole in 
the scattering matrix which can correspond to a bound 
state (2^. In particular, a very strong final state inter- 
action (FSI) is observed in the dd — > 4 He77 reaction close 
to kinematic threshold and is interpreted as a possible 
indication of 4 He — 77 bound state [23|. This suggests, 
that the 4 He — 77 system is a good candidate for the ex- 
perimental study of a possible binding. This conclusion 
is strengthened by the predictions in Reference [lj. 



The experiment was based on the measurement of the 
excitation function of the dd 3 Hep7r~ reaction for en- 
ergies in the vicinity of the 77 production threshold and on 
the selection of events with low 3 He center-of-mass (CM) 
momenta. In the case of existence of the 4 He — 77 bound 
state we expected to observe a resonance-like structure 
in the excitation function below the threshold for the 
production of the 4 He — 77 system. The structure would 
indicate a dd — > ( 4 He — r\)bound — > 3 Hep7r~ reaction ap- 
pearing on the continuous background originating from 
the direct dd — > 3 Hep7r~ process. 

We expect that the decay of such a state proceeds 
via the absorption of the 77 meson on one of the nucle- 
ons in the 4 He nucleus leading to the excitation of the 
N* (1535) resonance which subsequently decays in pion- 
nucleon pair. The remaining three nucleons play the role 
of spectators and they are likely to bind forming a 3 He or 
3 H nucleus. In the case of a similar system, the ^He hy- 
while ongo- pernucleus, it was observed that in the ir~ decay channel 



HHH, the theoret- 



However, as stated in Reference 
ical predictions for the width and binding energy of the 
77-mesic nuclei are strongly dependent on the subthresh- 
old 77-nucleon interaction which is not well understood. 
Therefore, direct measurements which could confirm the 
existence of the bound state are mandatory. 

Taking into account the above arguments and the fact 
that in the light nuclei systems the bound states are ex- 
pected to be much narrower compared to the case of the 
heavy nuclei [26^ . we performed a search for 77-mesic 4 He 
at the Cooler Synchrotron COSY-Julich with the WASA- 
at-COSY detector 1271. 



the decay mode ^He — > 3 Hep7r~ is dominant [28j. 

According to the discussed model, there exist four 
equivalent decay channels of the ( 4 He — rj) bou nd state: 



• ( 4 He - rj) bound -> 3 Hep7r" 



( 4 He - r)) bound -> 3 Hen7r° 



• ( 4 He - r)) bound -> 3 Hp7T° 



• ( 4 He - rf) bound -> 3 Hn7r H 



In the reported experiment we concentrated on the 
3 Hep7r~ decay mode. 

The WASA-at-COSY detector is described in detail in 
Ref. [29L l3o| . It consists of two main parts: the For- 
ward Detector dedicated to the measurement of forward- 
scattered projectiles and target-recoils, and the Central 
Detector, optimized for measuring of photons, electrons 
and pions originating from decays of mesons and excited 
baryonic states. The forward part consists of several lay- 
ers of plastic scintillators allowing for particle identifica- 
tion on the basis of the AE-E and AE-AE information 
and a proportional drift chamber providing track coordi- 
nates. The Central Detector is composed of an electro- 
magnetic calorimeter, a cylindrical drift chamber and a 
barrel of plastic scintillators. A superconducting solenoid 
provides a magnetic field for momentum determination 
of the tracks of charged particles measured in the drift 
chamber. The scintillators provide fast signals for the 
first level trigger, and together with the drift chamber 
and the calorimeter, are used for charged particle identi- 
fication via AE-p and AE-E methods. WASA-at-COSY 
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Edep(FRH2) [GeV] 



FIG. 1. Experimental spectrum of the energy losses in first 
two layers of the Forward Range Hodoscope. The area used 
for 3 He identification is indicated by the red line. The empty 
area below 0.05 GeV in the Edep(FRHl) distribution is due 
to the preselection condition. The regions corresponding to 
protons, 3 He and 4 He are clearly visible. 



uses an internal target system which provides pellets of 
frozen hydrogen or deuterium. During the present experi- 
ment the cooling system of the superconducting solenoid 
was broken and, therefore, no magnetic field was pro- 
vided. 

During the experimental run the momentum of the 
deuteron beam was varied continuously within each ac- 
celeration cycle from 2.185 GeV/c to 2.400 GeV/c, cross- 
ing the kinematic threshold for r\ production in the 
dd — > 4 He?/ reaction at 2.336 GeV/c. This range of beam 
momenta corresponds to a variation of 4 He — 77 excess 
energy from -51.4 MeV to 22 MeV. 

The identification of the 3 He was conducted using the 
AE — AE technique, comparing the energy losses in two 
layers of the Forward Range Hodoscope (Fig. [I}. 

The energy loss in the Plastic Scintillator Barrel was 
combined with the energy deposited in the Electromag- 
netic Calorimeter to identify protons and pions (Fig. [2]). 

The outgoing 3 He nucleus plays the role of a spec- 
tator and, therefore, we expect that its momentum in 
the CM frame is relatively low and can be approximated 
by the Fermi momentum distribution of nucleons inside 
the 4 He nucleus. This signature allows us to suppress 
background from reactions leading to the 3 Hep7r~ final 
state but proceeding without formation of the intermedi- 
ate ( 4 He — rebound state and, therefore, resulting on the 
average in much higher CM momenta of 3 He (see Fig. [3]). 

Therefore, we compare the excitation functions for the 
dd — > 3 Hep7r~ reaction from the "signal-rich" region cor- 
responding to the 3 He CM momenta below 0.3 GeV/c and 
the "signal-poor" region for the 3 He CM momenta above 
0.3GeV/c. The number of events as a function of the 
beam momentum is shown in Fig. 21 At this stage of 
the analysis the excitation function is not normalized to 
the luminosity and it is not corrected for reconstruction 
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FIG. 2. Comparison of the Monte Carlo simulation (upper 
plot) and the experimental spectrum (bottom plot) of the en- 
ergy loss in the Plastic Scintillator Barrel (x-axis) combined 
with the energy deposited in the Electromagnetic Calorime- 
ter (y-axis). The green and red curves represent the applied 
graphical condition to separate protons and pions. 



efficiency. The obtained functions for both regions are 
smooth and no clear signal, which could be interpreted 
as a resonance-like structure, is visible. We checked also 
for possible structures in the difference between the ex- 
citation functions for the "signal-rich" and "signal-poor" 
region. We multiplied the function for the "signal-poor" 
region by a factor chosen in such a way, that the differ- 
ence of the two functions for the second bin of Q is equal 
to zero. This difference is presented in Fig. @] (bottom) 
in order to examine the shape of the excitation function 
before any further selection criteria are applied. The ob- 
tained dependence is flat and is consistent with zero. No 
resonance structure is visible. 

In addition, further observables were taken into ac- 
count in order to reduce the background. Additional se- 
lection criteria on the p and 7r~ kinetic energy distribu- 
tions and the p — tt~ opening angle in the CM system 
were applied. In the N* rest frame this angle is exactly 
equal to 180° but due to the Fermi motion it is smeared 
by about 30° in the reaction CM system (see Fig. [5]). We 
also applied a condition on the relative p — ir~ angle in 
the CM system in the range of (140°-180°). 

The experimental spectra of CM kinetic energies 
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FIG. 3. (Upper plot) Distribution of the 3 He momentum in 
the CM system simulated for the processes leading to the cre- 
ation of the 4 Hery bound state: dd — > ( 4 Hen) bound — > 3 Hep7r~ 
(red area) and of the phase-space dd — > 3 Hep7r~ reaction 
(black line). The simulation was done for a momentum of 
the deuteron beam of 2.307 GeV/c. The Fermi momentum 
parametrization was taken from [3l| |. (Lower plot) Experi- 
mental distribution of the 3 He momentum in the CM system. 
In both plots the dashed line demarcates the "signal-poor" 
and the "signal-rich" regions. The decrease of the counts at 
0.48 GeV/c is due to the geometry of the border of the barrel 
and the end-caps of the Scintillator Barrel detector which was 
used in the p — n~ identification process. This region has no 
relevance in the next steps of the analysis. 



of protons and pions are compared in Fig. [S] to the 
distribution expected for the signal reaction dd — ► 
( 4 He — f])i, oun d — > 3 Hcp7r~. For further analysis we se- 
lected the kinetic energy of protons smaller than 200 MeV 
and of pions from the (180, 400) MeV interval. 

After the application of the described conditions the 
number of selected events in each excess energy (Q) in- 
terval was divided by the corresponding integrated lumi- 
nosity and corrected for the reconstruction efficiency. 

The absolute value of the integrated luminosity in the 
experiment was determined using the dd — > 3 Hen re- 
action and the relative normalization of points of the 
dd —± 3 Hep7r~ excitation function was based on the quasi- 
elastic proton-proton scattering [32] . 

The luminosity as a function of the excess energy is 
shown as triangles in Fig. is flat within the statisti- 
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FIG. 4. Excitation function for the dd — > 3 Hep7r _ reaction 
for the "signal-rich" region corresponding to 3 He momentum 
below 0.3GeV/c (upper panel) and the "signal-poor" region 
with 3 He momentum above 0.3GeV/c (middle panel). Dif- 
ference of the excitation functions for the "signal-rich" and 
"signal-poor" regions after the normalization to the second 
bin of Q is shown in the lower panel. The black solid line rep- 
resents a straight line fit. The threshold of 4 He — 77 is marked 
by the vertical dashed line. 



cal uncertainties. The geometrical acceptance is about 
60% and the overall efficiency including all selection con- 
ditions applied in the analysis is about 18% along the 
whole excess energy range. It is important to stress that 
both acceptance and efficiency are smooth and constant 
over the studied range. 

The excitation function obtained after the selection cri- 
teria on energy and opening angles, the correction for the 
efficiency and the normalization to the luminosity is pre- 
sented in Fig. [5J It can be well described by a second 
order polynomial (dashed line) resulting in a chi-squarcd 
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FIG. 5. (Upper panel) Simulated distribution of the p — n~ 
opening angle in the CM system for dd — > ( 4 Hej7)b olln[ j —¥ 
3 Hep7r~ reaction (red histogram) and for the phase-space 
dd —>■ 3 Hep7r~ reaction (blue line). (Lower panel) Experimen- 
tal distribution of thep — n~ opening angle in the CM system. 
In both plots the red dashed line separates the "signal-poor" 
and the "signal-rich" regions. 



FIG. 6. Kinetic energy distribution of protons (top plot) and 
7T~ (bottom plot) in the CM frame obtained from experiment 
(points) and from the simulations of a signal reaction (lines). 
The red dashed line indicates the boundary of the applied se- 
lection criteria: (Ekin° M < 200 MeV, Ekin° M G (180,400) 
MeV). Please note that the ranges on the X axes are different. 



value per degree of freedom of 0.98 and slightly worse by 
a straight line (solid line) . As in the intermediate stage of 
the analysis (Fig. [2| , the final excitation function exhibits 
no structure which could be interpreted as a resonance 
originating from the decay of the ?7-mesic 4 He. 
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UPPER LIMIT FOR THE 

( 4 Her;) 6ound ->■ 3 Hep7r- CROSS-SECTION 



Since no signal originating from the formation of the 
4 He — r\ bound state was observed, we estimate an upper 
limit for its production via the dd — > (^Hei])t, oun( i — > 
3 Hep7r~ reaction. 

We assumed that a signal from the bound state in the 
excitation curve determined as a function of the excess 
energy Q with respect to the 4 He — r\ threshold can be 
described by a Breit-Wigner shape: 



<t(Q,E BE ,T,A) 



(i) 



{Q-E BE y + {^ 

where Ebe is the binding energy, T is the width of the 



bound state and A is the cross-section at the central en- 
ergy (Q = E B e) for the dd -> ( 4 He - rj) bound -> 3 Hep7r" 
reaction. In this way, we assume that there is no inter- 
ference between the signal and the non-resonant back- 
ground. In order to determine an upper limit for the 
cross-section for formation of the He — rj bound state 
and its decay into the 3 Hep7r~ channel we fitted the ex- 
citation function with a polynomial describing the back- 
ground (first and second order) combined with the Breit- 
Wigner function. In the fit, the polynomial coefficients 
and the amplitude A of the Breit-Wigner distribution 
were treated as free parameters. The binding energy Ebe 
and the width T were fixed during the fit. 

The fit was performed for various values of the binding 
energy and the width representing different hypothesis 
of the bound state properties. The binding energy and 
the width were varied in the range from to -30 MeV 
and from 5 to 35 MeV respectively. In each case, the 
extracted value of A is consistent with zero within the 
statistical uncertainties, which confirms the hypothesis 
of non-observation of the signal. 

In order to calculate an upper limit for the dd — > 
( 4 He?7) bound 3 Hep7r~ cross-section, the standard de- 
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FIG. 7. Geometrical acceptance (full black squares), overall 
efficiency (open red circles) and luminosity (full blue trian- 
gles) as a function of the excess energy. The right axis denotes 
the luminosity. 
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TABLE I. The upper limit for the cross-section for the 
bound state formation and decay in the process dd — >• 
( 4 He — r))b oun d —¥ 3 Hep7r~, determined on the 90 % confi- 
dence level. The values were obtained from a fit of a Breit- 
Wigner function combined with first and second order polyno- 
mials (aiin and a qU ad respectively) with different fixed values 
of binding energy Ebe and width F. 
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FIG. 8. Experimental excitation function for the dd — > 
3 Hep7r - reaction obtained after the normalization of the 
events selected in individual excess energy intervals by the 
corresponding integrated luminosities. The dotted and solid 
lines correspond to the second and the first order polynomials 
fitted to the data. 
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FIG. 9. Upper limit at 90 % confidence level of the cross- 
section for formation of the 4 He — r\ bound state and its decay 
via the dd — > ( 4 He — rj)bound — > 3 Hep7r~ reaction as a function 
of the width of the bound state. The binding energy was set 
to Ebe=-2Q MeV. The grey area at the bottom represents the 
systematic uncertainties. 



SYSTEMATICS 



viation of the A values (<ta) obtained from the above 
described fit were multiplied by the statistical factor k 
equal to 1.28155 corresponding to the probability confi- 
dence level (CL) of 90%. The final results were obtained 
by averaging the upper limits derived from fits with a 
background described by the second and first order poly- 
nomials. 

The examples of the obtained upper limits are given 
in the last column of Table I. One can notice that these 
limits depend mainly on the width of the bound state 
but only slightly on the binding energy. The result for 
E BE =-20 MeV is shown in Fig. EE 



Systematic checks were performed by studying the sen- 
sitivity of the result to the variation of the selection con- 
ditions performed in the analysis and the assumption 
taken in the fitting procedure. 

Changing the range of the above described selection 
conditions within ± 10 % gives a result consistent within 
the statistical uncertainties. 

The smooth reconstruction efficiency and the luminos- 
ity dependency as a function of excess energy is of high 
importance because it eliminates the possibility of the 
creation of an artificial signal due to fluctuation of the 
acceptance or the luminosity. 

Two methods were applied to extract the luminosity 
dependency as a function of Q. In addition to the nor- 
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malization calculated on a bin by bin basis, we have es- 
timated the luminosity dependence of Q using a fit of a 
first order polynomial to the data. The results of both 
methods are in agreement. However, an overall normal- 
ization uncertainty of luminosity is equal to 11.5 % [32| 
and this value is one of the contributions to the system- 
atic uncertainty of the upper limit. 

The description of the background shape with 
quadratic and linear functions produces additional sys- 
tematic uncertainty, which was estimated as: 

^ ip quad &lin) 

2 

The systematic error grows almost linearly with the as- 
sumed bound state width from about 5 % (r = 5 MeV, 
E BE = -20 MeV) to 33 % (T = 35 MeV, E BE = -20 MeV) 
and we take that range as an estimate of the systematic 
uncertainty due to the assumed shape of the background. 

An important source of systematic errors comes from 
the Fermi momentum distribution of nucleons inside the 
4 He nucleus applied in the simulations. We adapted 
the Fermi momentum parametrization described in [3l|, 
which is derived from the work of McCarthy, Sick and 
Whitney (33|. However, as it is shown in [34| . the alter- 
native derivation of the Fermi momentum distribution 
done by Nogga [35] is narrower by about 25 %. Even if 
the choice of the given Fermi model does not influence the 
experimental method, it affects the overall acceptance of 
the 3 He ions in the Forward Detector, and adds an addi- 
tional systematic error of 8 %. 

In principle, a complete description would require the 
application of a momentum distribution with the embed- 
ded N* resonance. However, up to now, such a descrip- 
tion of the momentum distribution is unavailable. There- 
fore, we approximate this distribution by Fermi momen- 
tum distribution of nucleons inside the 4 He nucleus. 

Adding the above estimated contributions in quadra- 
ture we obtain systematic uncertainty of the upper limit 
which grows with bound state width from 15 % to 36 %. 



SUMMARY 

We performed a search for the 4 He — rj bound state 
via exclusive measurement of the excitation function for 
the dd — > 3 Hep7r~ reaction. The measurement was car- 
ried out with the internal deuteron beam of the COSY 
accelerator scattered on a deuteron pellet target and 
with the WASA-at-COSY detection system used for reg- 
istration of the reaction products. During the experi- 
mental run the momentum of the deuteron beam was 
varied continuously within each acceleration cycle from 
2.185 GeV/c to 2.400 GeV/c, crossing the kinematic 
threshold for r\ meson production in the dd — > 4 He rj re- 
action at 2.336 GeV/c. This range of beam momenta 



corresponds to an interval of the excess energy in the 
4 He - r] system from -51.4 MeV to 22 MeV. 

For the first time in the experimental search for mesic 
nuclei all ejectiles were measured and the reaction was 
identified exclusively. 

No signal from 77-mesic 4 He was observed. The up- 
per limit for the cross-section for the bound state forma- 
tion and decay in the process dd — > ( 4 He — rj)i, oun d — >• 
3 Hep7r~, was determined on the 90 % confidence level 
and it varies from 20 nb to 27 nb for the bound state 
width ranging from 5 MeV to 35 MeV, respectively. The 
upper limits depend mainly on the width of the bound 
state and only slightly on the binding energy. 
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